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Abstract 

The break-up of the deuteron during deeply-virtual Compton scattering, j*d — ► j^np, is explored. 
In the effective field theory describing nucleon dynamics at momenta below the pion mass, the EMC 
effect results from four-nucleon interactions with the twist- 2 operators, appropriate for describing 
forward, and near-forward, matrix elements in the two-nucleon system. We point out that the 
break-up of the deuteron to low-energy final states during deeply-virtual Compton scattering is 
a process with which to explore strong-interaction physics closely related to that responsible for 
the EMC effect. The single-nucleon contribution to the break-up depends on the moments of the 
spin-dependent structure functions and contributions from local four-nucleon operators. Experi- 
mental deviations from the single-nucleon prediction would provide a probe of strong interactions 
complimentary to the EMC effect. 
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A naive glance at QCD would suggest that the nucleus is a non-perturbative collection 
of quarks and gluons with no discernible structure. However, this is not the case, and a 
remarkable aspect of nuclear physics is that nuclei may be viewed as a collection of relatively 
weakly-interacting nucleons. The small binding energy of the deuteron implies that it very- 
nearly resembles a free proton and neutron; its magnetic moment is very close to the sum of 
proton and neutron magnetic moments, and its structure functions are very close to the sum 
of the nucleon structure functions. Deviations from the single-nucleon structure functions, 
known as the EMC effect 1 , were experimentally observed 0] in nuclei some twenty years ago 
and caused a great deal of excitement, e.g. 0,0,0,0- However, as the nucleons do interact, 
albeit weakly, one expects matrix elements of all electro-weak probes to be modified from 
their single-nucleon values, and the EMC effect should be considered to be nothing more than 
another demonstration of this. Significant efforts have gone into quantitatively explaining 
the EMC effect in explicit models of nuclei, e.g. Ref. p,l3]. One would like to uncover the 
dominant strong-interaction mechanism that allows for a quantitative explanation of the 
effect, but given that one cannot accurately predict the deuteron quadrupole moment from 
modern potential models 2 to better than ~ 5%, an approximate explanation would suffice 
at present. 

Within the framework of the pionless effective field theory, EFT(^), it is easy to under- 
stand and describe the EMC effect. By constructing all the operators required to reproduce 
the matrix elements of the twist- 2 operators in multi-nucleon systems, one sees that oper- 
ators involving more than one nucleon are not forbidden by the symmetries of the strong 
interaction, and therefore must be present. While observation of the EMC effect twenty 
years ago may have been surprising to some, in fact, its absence would have been far more 
surprising. In this work we write down the most general set of operators describing the 
interaction between two nucleons in an S-wave that reproduce the matrix elements of the 
twist-2 operators. As the origin of the EMC effect is strong-interaction physics at distance 
scales of order the Compton wavelength of the pion or smaller, there will be an analogous 
deviation from the single-nucleon prediction for transitions between two nucleons in the 3 5*i 
channel and the ^So channel. Therefore, we suggest that an experimental determination 
of the break-up of the deuteron during Deeply- Virtual Compton Scattering (DVCS), with 
very low-energy final-state nucleons, provides a probe of a purely-nuclear effect which is 
complimentary to the EMC effect. 

In any given hadron, the moments of the parton distribution functions are given by 

I dxx n - x (q(x) + (-l) n q(x)) ; 
Jo 

dxx" 1 - 1 (Aq(x) + Aq{x)) , (1) 

where q (q) is the quark (antiquark) spin-averaged distribution and Aq (Aq) is the helicity 



1 The expression "EMC effect" usuaily refers to deviations from single-nucleon structure in nuclei with 
A ^S> 2. Here we apply the expression to the deuteron and therefore, in the present context, the EMC 
effect is relevant to extractions of the neutron structure functions. 

2 The quadrupole moment of the deuteron is successfully described in effective field theory by the inclusion 
of a local operator at next-to-leading order in the expansion Q . The discrepancy between potential model 
calculations and the experimental value is due to the omission of physics at distance scales smaller than the 
Compton wavelength of the pion, which may or may not be describable by a meson-exchange mechanism. 
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distribution (we do not consider the transversity distribution). From QCD, these moments 
have the clear interpretation as the forward matrix elements of local operators, 



{ vX,;n = (^r 1 ?^ D, 2 ... D, n} q 

= ®r X ll{K D, 2 ... D, n} r b q 

O^'l-vn = (iT^Ql^ D» 2 ... D, n} 75 q 

0^i.., n = {^ll^ D, 2 ... D, n} J 5 T b q , (2) 

where the indices enclosed by {...} are symmetrized and their traces are removed. An 
analysis of matrix elements of these operators in the single-nucleon sector has been carried 
out at the one- loop level in the chiral expansion 0, 11, 0, k|. Those works have 



been extended from QCD to quenched QCD [3| and partially-quenched QCD [lil, 1(| for 
applications to lattice QCD simulations. In addition to forward matrix elements of these 
operators, off-forward matrix elements have been considered |17| . and the issue of soft-pion 
production during DVCS has been explored [3, [l^ 3 . For very low-energy processes, those 
involving momentum and energy transfers much less than the pion mass, one can work in 
EFT(jf). In this energy regime, the matrix elements of the operators in eq. (J2J) between 
single-nucleon states are reproduced by 

{ v^ - Mn ~ l (^"V - *W NN + ... 

^ } L»n - Mn ~ l i^' 1 )^ v^Vh - tw Nr b N + ... 

Ofi.., n - M n ~ l (x"- 1 )^) v {fll v, 2 ... v, n _, NS, n} N + ... 

Ofi^ - M n ~ l (ar-\<p v { ^ 2 ... v, n _, NS, n} r b N + ... , (3) 

where (x p ) (oj and (x p ) (o) are the p'th moments of the isoscalar vector and axial-vector 
PDF's in the nucleon, respectively, and (x p )^i) and (x p ) are the p'th moments of the 
isovector PDF's in the nucleon, respectively. M is the nucleon mass. For n = 1 baryon- 
number conservation gives (l) g (o) = 3; isospin-conservation gives (l) 9 (i) = 1; and matching 

to the axial-vector currents gives (1) (oj = go, the singlet axial charge 4 , and (1) (i) = gA, the 

q A q A 

isovector axial charge. The nucleon four-velocity is v , and in the nucleon rest frame becomes 
v = (1, 0, 0, 0). The ellipses denote terms that are higher order in the derivative expansion, 
suppressed by powers of p/m n . 

When considering the parton distributions of the deuteron, the simplest nucleus, inser- 
tions of the operators of eq. (j3J) defines the impulse approximation. Given the operator 
structure, it is clear that one simply recovers the naive sum of contributions from the proton 



3 The errors made in Ref. [18j and addressed in Ref. |l9j have been replicated in Ref. [20( . This latter paper 
does not properly account for the fact that chiral symmetry is spontaneously broken, and thus leads to 
results — like gA = 1 — which, while possibly interesting, are not relevant to QCD. 

4 The "q" subscripts on the operator coefficients in eq. 10) and eq. (@J indicate that we are considering matrix 
elements of the quark operators only, and not the purely-gluonic operators -including the anomaly- which 
mix with the quark operators under renormalization. The singlet axial charge defined here, g$, is therefore 
scale dependent. 
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and neutron. However, it is also clear that, in the effective field theory expansion, there are 
contributions from higher-dimensional operators involving more nucleons. The leading-order 
interactions involving two-nucleons are 



e 



(n),0 
V, Hx-Hn 



M 



n-l 



(n) 
q 



a(n),a 

U V, fil-Hn 

n(n),0 

U A, (Mi. .fl n 

n(n),a 
°A, Hl-Hn 



... Vfln} (iV T PfiV) t N T P?N 
+ $} v { ,iv, 2 ... v, n} (N T P«Nf N T P«N 
+ 0$} ... v, n (N T PtN) j N T Pt } N 

(n t p£nV n t p?n 



M n-l 7 (n) l£ abc 



Vn} 



M n ~ l 



V {»l V H2 



7 , 7 >Mn}a/37 ? , 

U IX n -l 6t U - 



M n ~ l v {l „ r 



Ml 



(N T PiNy N T PJN 
1 (N T P?Ny N T Pt } N + h.c. , 



where with n 



1 we have /3<£g 



(4) 



conservation gives 
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P0;q 



the axial-vector current gives p 



(i) 
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(due to available Lorentz indices) and baryon-number 

'' The objects Pf and 



0; isospin-conservation gives 7^ = 0; and matching to 



Li, 



21 



-2L 2A and 

P3 [221 l2ii |24| are the spin-isospin projectors for the and 3 Si channels, respectively, 



pa 



V2 



r 2 r 



So 



pM 
*3 



V2t 2 ®S 2 S» 



(5) 



where is the usual spin operator which acts on the nucleon fields. The contributions 
of such operators to the parton distributions of the deuteron can be computed straightfor- 
wardly in EFT(^f) 0, 6 . However, it is convenient to switch notation to that of EFT(^) 
constructed in terms of di-baryon fields [2U dEFT(^), in which the leading effects of 
the finite effective range are resummed. 

In this work we will discuss relevant aspects of the di-baryon formalism, dEFT(^) [2f 



271 . l28| . In terms of nucleon and di-baryon degrees of freedom, the leading-order low-energy 



strong interactions for |p| <C rn n /2 are described by a Lagrange density of the form 

V 2 



c 



+2/3 



id + 



2M 



N + t 



V 2 

ido + TTT - ^3 



f - S 



t\ N T P3N + h.c. 



Vi 



AM 

si N T P?N + h.c. 



V 2 



(6) 

where N is the nucleon annihilation operator, we have used the four-vector = (0,tj) with 
P the 3 Si di-baryon annihilation operator with spin- index j, and s a is the l So di-baryon 
annihilation operator with isospin- index a. In the Lagrange density in eq. (jUJ), a factor of 



5 The cross-sections for the weak disintegration of the deuteron and other two-nucleon weak processes 
depend on the parameter Li t A- The axial operator with coefficient £2,^4 contributes — together with the 
sum of the neutron and proton contributions — to the strangeness content of the deuteron. 

6 The effective field theory can be constructed in terms of operators that do not vanish on the particle 
mass-shells. This is true for any consistent theory, and operators that only contribute off mass-shell can 
always be eliminated. 
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FIG. 1: Diagrams in dEFT(^) that provide the leading-order contribution to the matrix elements 
of the twist- 2 operators in the deuteron. The thick solid lines with solid circles at their midpoint 
correspond to fully- dressed di-baryons, while the thin lines correspond to nucleons. The crossed 
circles correspond to insertions of the twist- 2 operators. 



the nucleoli mass M (2M) has been absorbed into the definition of the nucleon fields (di- 
baryon fields). The S-wave interactions are enhanced by a factor of the expansion parameter, 
1/Q, and are treated non-perturbatively. However the interactions that induce mixing with 
higher partial waves, e.g. 3 S% — 3 Di mixing, are suppressed by at least Q 2 and so only S-wave 
to S-wave interactions are required to the order we are working. To recover the scatte ring 
amplitudes in both S-wave channels, the constants that appear in eq. © are found to be [26j] 

2 8vr 2 8vr a 2 f 1 \ 2/1 \ _ 

where /i is the renormalization scale, 03 and r% are the scattering length and effective range 
in the 3 Si channel, and a\ and r\ are the scattering length and effective range in the ^So 
channel. Insertions of the twist-2 operators give rise to operators of the form 



+ Po£ ^1^2 ••• *w tie + p$ u {WJ v w ... v^ n 



a V, flL.fln ^ iV1 % l£ V {»1 V H1 ■■■ V Hn} S b S C 

- M- 1 v {fll v, 2 ... ie**** v a t\ t, 

- M- 1 ?W v {fll v, 2 ... [ 4 + h.c. ] . (8) 

With n = 1 we have (3^1 = and baryon-number conservation gives fiW = — /3q^ = 6; 
isospin-conservation gives = 2; and matching to the axial-vector current gives = 
— 2^2,A/2M v /rir 3 and ajp = £i i A/2M v /F7r 3 " 0. (Recent work has shown that it is possible 
to compute the coefficients of the operators in eq. (JHJ) in lattice QCD by making use of 
background-field methods j29l l3c|.) 

It is now straightforward to determine the matrix elements of the relevant operators 
between deuteron states. Keeping only the AL = contribution from the diagrams shown 
in Fig. ^ the matrix elements of the isoscalar- vector quark operators at momentum transfer 
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q <C m n are 



(d\ { y>Z., n \d) = M n ~ l v {fll v, 2 ... v, n} 



2(x- 1 ) <? (o ) F c (|q|) - < r 7?:; 



7^3 

where the expression is given in terms of the electric charge form factor 

7^3 



(9) 



*b(|q|) 



7 r 3 



q r 3 




tan" 1 I — I - 1 



(10) 



The single-nucleon contribution has been made explicit by writing 



= 2(x"- 1 ) 9 (o ) + < . (11) 

Here we see explicitly that the contribution from the nucleon structure functions (first term 
in the square brackets of eq. ©) is accompanied by a nuclear contribution (second term in 
the square brackets of eq. (jOJ)), making clear the origin of the EMC effect. Notice that for the 
baryon-number current (n = 1), fj^l = and the matrix element of eq. © is entirely given 
by the leading-order electric charge form factor of the deuteron; as emphasized in Ref. |28l ]. 
it is a peculiar feature of dEFT(^) (unlike EFT(^)) that higher-dimensional local di-baryon 
operators contain contributions that should be interpreted as leading-order in the effective 
field theory expansion. Several points are worth emphasizing as they bear on conventional 
explanations of the EMC effect: 1) relativistic corrections appear at the next order in the 
effective field theory expansion and can be systematically included. These corrections have 
been computed for various deuteron properties in EFT(^) and found to be small [1.8]; for 
static processes these corrections scale as 7 2 /M 2 . Notice that this expansion parameter 
is suppressed relative to the expansion parameter 7/777.^ characteristic of the size of the 
purely nuclear effects in EFT(^f) and dEFT(^f). 2) "Fermi motion" (or more appropriately, 
"zero-point energy") of the deuteron, is included in EFT(/) and dEFT(^); terms of the form 
p 2 /2m where p is a nucleon momentum are summed to all orders through the non-relativistic 
nucleon propagators from eq. (JSj). 3) As pointed out above, in EFT(^) and dEFT(^f) D-state 
admixtures enter at higher order in the effective field theory expansion (if. 

There is also a "tensor EMC effect". This is induced by the operator in eq. (jHJ) with 
coefficient /^-g , and corresponds to transferring AL = 2 to the hadronic system. This does 
not contribute to DIS, as it is not a forward matrix element, but will contribute to DVCS. 
If the deuteron polarization in the final state could be measured after DVCS on a polarized 
deuteron target, then in principle this object could be measured. Experimentally, such a 
measurement will not be easy. 

A different situation in which the effects of the interactions between nucleons during 
DVCS will lead to modifications from the impulse approximation is the break-up of the 
deuteron 7 . Such deviations are, in some sense, intimately related to the EMC effect as they 
result from different four-nucleon operators connecting S-wave states. They are analogous to 
the local four-nucleon operators which cause deviations from the impulse approximation for 



7 For a recent discussion of the process e + d — > e + n + p in the 1.5 < Q 2 < 4 GeV 2 range, see Ref. 1311. 
Deviations from the impulse approximation in DVCS on nuclei has been considered recently in Refs. [32( 
and 0. 
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FIG. 2: Diagrams in dEFT(^) that provide the leading-order contribution to the matrix elements 
for the break-up of the deuteron during DVCS. The thick solid lines with solid circles at their 
midpoint correspond to fully-dressed di-baryons, while the thin lines correspond to nucleons. The 
crossed circles correspond to insertions of the twist- 2 operators. 



low-energy np — > d'-f (which result from interactions between nucleons that are not related 
by gauge- invariance to the nucleon-nucleon scattering amplitude). It may or may not be the 
case that pion-exchange currents can reproduce much of this contribution. 

It is straightforward to show that for very-low energy and momentum transfers to the 
hadronic system, the amplitude for DVCS-dissociation of the deuteron can be written as 



A 



X M 



n—l 



(d) 

V Vn-l £ fi n } 



(12) 



for the neutron and proton in a final-state, where is the deuteron polarization vector, 
and Un tP are the spinors associated with the neutron and proton respectively. For the sake of 
simplicity, we give an expression for the coefficient X in eq. ()12j) . resulting from the diagrams 
in Fig. |2l for kinematics where the energy-transfer and momentum-transfer to the hadronic 
system are equal, i.e. light-like. We find 



X 



2i 4 



1 1 


7 3 y/1 - 7r 3 




2 -i 
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/ 



n-1 
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P + 7 
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1 1 , 
ai 2 



(13) 
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where p is the momentum of one of the nucleons in the final state 8 . While the spin- 
dependent structure functions determine the single-nucleon contributions to the deuteron 
break-up (first term in the square brackets of eq. (|T3j) proportional to (x n ~ 1 ) w ), there is 

an additional contribution coming from short-distance strong interactions (second term in 
the square brackets of eq. (|13|) proportional to c^ Tl * ) ). This amplitude is similar in form to 
that of np — > j34l 35j ] for obvious reasons 
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We emphasize that the validity of dEFT(^f) requires that the energy imparted to the 
hadronic system be significantly less than m n . For higher energies, pions must be included 
explicitly in the description. Even with the inclusion of pions as dynamical fields, local 
four-nucleon operators will continue to contribute to nuclear matrix elements of the twist- 2 
operators. In the pionful theory these operators describe contributions from physics at scales 
of the order of and smaller than the chiral symmetry breaking scale. It is important to point 
out that the same two-body operators that are inserted in the deuteron will contribute to 
DVCS on the triton and 3 He, with three-body operators expected to be suppressed by their 
naive engineering dimensions |36L |3?| . 

In this work we have discussed the EMC effect in terms of the matrix elements of local 
operators. We have made no attempt to determine the underlying mechanism; rather we 
have demonstrated how to describe the effect in a model-independent manner. As for the case 
of the single nucleon, only positive moments of the parton distributions can be related to the 
forward matrix elements of local operators, and so it is natural to consider only the positive 
moments of the parton distributions of nuclei, including the deuteron. By contrast, the small- 
x behavior is not related to forward matrix elements of local operators, and arises from long- 
distance strong-interaction contributions to the electromagnetic-nucleus interactions. So 
while the positive moments can be considered properties of nucleons and nuclei, the small-x 
behavior is a property of the interaction matrix element, and in general, is different for each 
probe. For forward processes and those involving only small momentum transfers we have 
used EFT(^) and dEFT(^) to compute the EMC effect and the DVCS-dissociation of the 
deuteron in terms of single-nucleon contributions and four-nucleon interactions due to short- 
distance physics. We consider this paper to be the first step toward a detailed calculation 
of the dissociation process with realistic kinematics. We estimate that a measurement of 
the break-up cross sections will show deviations from the single-particle estimate at the 
level of 5% — > 10%, in analogy with the process np — > dj, and in keeping with the EMC 
effect. Measurement of this break-up cross section at Jefferson Laboratory or elsewhere, and 
a determination of the deviation from the single-nucleon effect, would provide a fascinating 
probe of the strong interactions. 



8 Notice that while this form- factor sums all orders in |p|/7, it is perturbative in the expansion parameter 
\p\/m n . 

9 Both terms in the square brackets of eq. I|13|l naively appear to be of the same order in the expansion, 

- (n) 

however one power of the expansion parameter, Q, has been hidden in & q for ease of notation. 
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